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Abstract 

We suggest a simple model to describe the reversible field-induced transfer 
of a single Xe-atom in a scanning tunneling microscope, — the Eigler-switch. 
The inelasticly tunneling electrons give rise to fluctuating forces on and damp- 
ing of the Xe-atom resulting in an effective current dependent temperature. 
The rate of transfer is controlled by the well-known Arrhenius law with this 
effective temperature. The directionality of atom transfer is discussed, and 
the importance of use of non-equlibrium-formalism for the electronic environ- 
ment is emphasized. The theory constitutes a formal derivation and general- 
ization of the so-called Desorption Induced by Multiple Electron Transitions 

(DIMET) point of view. 
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Eigler and Schweizer [|l| have been able to manipulate Xe atoms and place them with 
atomic precision on a Ni surface using an ultra high vacuum scanning tunneling microscope 
(STM) operated at 4 Kelvin. Experiments by Eigler, Lutz and Rudge show [Q, that a single 
Xe atom physisorbed to a particular kink site on a single crystal Ni-(llO) surface, can be 
transferred reversibly between this surface-site and the W tip in the STM. These experi- 
ments are done at 4 Kelvin. Using voltage pulses of ±0.8V^ for 64msec, they are able to 
toggle the Xe atom from surface to tip and back to the same position on the surface at will, 
and find that the direction of transfer is the same as that of the tunneling electrons. For 
a particular tip-surface configuration with the Xe adatom on the surface, the transfer rate, 
r^^, goes as j4-9±o.2^ j ];3gij2g the tunneling current. The voltage range in these measurements 
is from 18mV to ISOmV with V/I = 906kQ ± 2%. These phenomena has been investigated 
theoretically by Walkup, Newns and Avouris and Gao, Persson and Lundquist 0. These 
authors suggest that the main mechanism behind the transfer is that the current excite the 
Xe atom vibrationally in the double-well potential, sustained by the van der Waal attraction 
to surface and tip. On the other hand the Xe atom dissipates energy to the surface phonons, 
so they calculate the rate of transfer using rate-equations (Pauli master equations) including 
the competition between dissipation to surface-phonons and "heating" by inelasticly tun- 
neling electrons. In Ref. P] the possible mehanisms behind the directionality of transfer is 
discussed, and they conclude that the adsorbtion- induced dipole is the dominant effect. 

In this paper we concentrate on showing that the key-features of the experiments can be 
explained from a less phenomenological point of view, based on a simple model. 

The main idea in our calculation is the following. We view the Xe atom as a quantum 
Brownian particle interacting with the environment of electrons in the tip and surface via 
the adsorbate-resonance, through which the electrons tunnel. Besides this environment 
the Xe atom also interact with an environment of surface-phonons. The influence of these 
environments on the atom is described in a path-integral framework using the influence- 
functional introduced by Feynman and Vernon |]^,|^, giving an effective action describing the 
motion of the atom. For two independent environments, the influence-functional will simply 
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be the product of the individual influence- funtionals [0. The influence- functional for a 
harmonic oscillator with a linear coupling to a continuous distribution of harmonic oscillators 
in thermal equilibrium has been calculated by Caldeira and Legget [0]. This influence- 
functional contains the fluctuating force and dissipation caused by the environment which in 
the classical limit are the ingredients in the Langevin equation describing the classical motion 
P|. Our contribution is the calculation of the influence-functional corresponding to the non- 
equilibrium electronic environment. In general, when a Brownian particle is in equilibrium 
with a heat-bath, the fluctuating force acting on the particle and the corresponding friction 
force are "connected" by the fluctuation-dissipation theorem (FDT). This is the situation 
when no chemical-potential difference between tip and surface drives a current through 
the adsorbate-resonanse, and the Xe atom is in equilibrium, the surface-phonons and the 
electrons acting as heat baths. When an external field in some way is transferring energy to 
the Brownian particle the relation between the fluctuating force and the dissipation, found in 
the equilibrium case, no longer holds. There will in general occur a extra current- dependent 
fluctuating force not "compensated" by friction. 

We show this general feature in our model, where the interaction with the inelasticly 
tunneling electrons give rise to the uncompensated fluctuating force. At sufficiently low 
temperature compared with the oscillation frequency in the adsorbtion well, as in the actual 
experiment, this force will be the dominant. The temperature independent dissipation is 
due to exitation of surface-phonons and electron-hole-pairs in the surface and tip by the 
vibrating Xe atom. The time scale of desorbtion is much longer than than the relaxation 
time of the dissipative system, thus we have a quasi-stationary situation a relaxation time 
after the current is turned on. In this quasi-stationary situation we can consider the Xe-atom 
as a damped harmonic oscillator with a wavepacket width given by the usual FDT-result 
plus a contribution due to the non-equilibrium fluctuating force. From this width we define 
our effective temperature. The desorbtion rate is then given by the well-known results for 
escape of a Brownian particle from a potential well in equilibrium with a heat-bath where 
the new effective temperature depends on the voltage. This leads to the experimentally 
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found result if we assume that the adsorbtion well contains about five bound states, as was 
done in Refs. 0] and . The dependence of the direction of transfer on the polarity of field 
follows in our model from the mean occupation of the Xe atom. We stress that it is crucial 
to use a non-equilibrium formalism to calculate the occupation. 

Let us now introduce the model and sketch the derivation of our results. The reason why 
a rare-gas atom like Xe actually can be imaged by STM is explained by Eigler et al. |0. Thus 
we will only consider this adsorbate state. Denoting the operator and energy associated with 
the Xe 6s-orbital d and ea, and the operators asociated with surface, respectively tip, Cfc, q, 
we write the Hamiltonian for the electronic part, 

Hei = Hsurf + Htip + ead?d + d\TkCk + T/q) + h.c. (1) 

k,l 

We define the weighted density of states for surface and tip, 

Ps{E)=Y. \n\H{E-eu) (2) 

k 

pt{E)=Y: mmE-e,). (3) 

For simplicity we use the wideband limit and take these to be constant. The current, 
corresponding to the model above, can be expressed by a non-equilibrium Green's function 
which is evaluated using the Keldysh technique [0, yielding 

m = (4) 

4e r°° duj {7ipt){Tcps){np{uj — eV) — np{uj)) 

h J-oo 271 {Tl{pt + Ps)y + {UJ - 

We have neglected the small influence of the vibrations of the Xe atom on the current. 

The parameters of the electronic Hamiltonian (|l]) depends on the position of the adatom. 
This of course leads to the picture of effective potential energy-surfaces, which is the basis of 
most surface dynamics calculations. We shall see that this dependence also is responsible for 
the fluctuations in the force experienced by the atom. Denoting by x the operator associated 
with the position of the Xe-atom, relative to the equilibrium distance from the surface, and 
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expanding the the onsite energy to first order in x we obtain the couphng between electrons 
and Xe 



Vel-vib = gid^d - {d^d)eq) X, (5) 

where g is the positive couphng-constant. (For simphcity we neglect the x-dependence of 
and Ti). Note that we have subtracted the equilibrium {eV = 0) mean population of the site, 
because this contribution is assumed already to be incorporated in the adsorbtion potential- 
surface. The physical picture of the interaction is that the temporary increase in charge 
on the site, when a electron tunnel through it, will cause the atom to be attracted more 
strongly towards the strongest image, which will be the one in the surface when the atom 
is located at the surface. Thus g will be positive. The interaction results in a fluctuating 
force acting on the Xe atom, due to the electrons, jumping on and off the site. This physical 
effect is what has been dubbed DIMET (desorbtion induced by multiple electron transitions) 



Tl| . On the other hand, the vibration of the atom will, because of this interaction, excite 
electron-hole pairs in the metal-surface and dissipate energy to the electronic part of the 
system. A similar model has been used by Persson and Baratoff to describe the influence 



of inelastic tunneling on the current in the STM. 

We use the Feynman- Vernon theory of infiuence-functionals to do a systematic calcula- 
tion of these features. Considering a quantum system interacting with a general environment, 
one is often only interested in the influence of the environment on the system. This can 
be accomplished by the reduced density operator, tracing out all superfluous information 
of the environment. This is elegantly done in a path- integral framework 0. The result of 
this procedure is an effective action of the system. This approach has been used on the 
problem of a quantum oscillator in equilibrium with a heat-bath consisting of a continuum 
of quantum oscillators, and is treated thoroughly in Refs. In our problem we will be 

concerned only with the behavior of the Xe-oscillator and thus trace out the equilibrium 
phonon environment as well as the non-equilibrium electronic environment. 

The propagator for the reduced density-operator for the oscillator can be written as a 



double path integral, 



J{xf,yf,t;xi,yi,0) = (6) 



where S is the action for the non-interacting oscillator, and J-'ei and J-'ph are the infiuence- 
functionals containing all relevant information of the electronic and phonon environments, 
respectively. These will, in general, couple the paths moving "forward" (x) and "backward" 
(y) in time. The phonon influence functional, corresponding to the oscillator coupled to a 
continuum of oscillators, has been evaluated in Ref. [0], 

^S^effix{T),y{T)) = -in \n[J^ph{x,y)] 

^ {x + y){x- y)dT (7) 
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The real part of Seff describes the friction [13] and the imaginary is the correlation function 



of the time-non-local fluctuating force. These two terms are related through the equilibrium 
FDT. 

The new feature here is the contribution to Sgjj from the non-equilibrium electronic 
environment. This is calculated using a technique similar to the one used in Ref. 
extended to cover this non-equilibrium case. The first order term is, 

ASfJ^'\x,y)= (8) 

- g {N{eV) - N{0)) f dt' {x{t') - 

Jo 

where N{eV^ is the mean population of the Xe site after the voltage is applied, and is given 
by, 

N[eV) = (9) 
°° Pt npjuj - eV) + ps upjuj) 
-oo {u - + TT^ipt + Ps^ 
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where the e^p^'^ = and e^^ = eV. This corresponds to a simple change of the potential 
energy surface due to the change of population of the site. The mean-population of the 
resonance increases for eV positive, and thus increase the energy-barrier between the surface 
and tip adsorbtion-sites, whenever the electrons tunnel from tip to surface. Reversing the 
polarity, so the electrons tunnel from surface to tip, decreases the barrier. When the Xe- 
atom is located on the tip, the roles played by tip and surface are exchanged. This can 
explain the observed dependence of transfer on the polarity. 

To second order in the coupling constant g we find terms similar to (|^) with a friction- 
coefficient, riei{eV), which can be expressed in terms of g, ps, pt and eV |ll5|. Besides these 



terms we find the additional imaginary term in ASl'jj'^\ 

{x{t) - y{T)) cos(cj(r - s))(a;(s) - y{s)) dr ds — , (10) 

with 

(^2 ^eV PsPt 
a(eV) = / — tt; rTTTT duj. (11) 

The effective action obtained so far is quadratic in the Xe-atom coordinate, and the 
remaining path-integral is Gaussian and can be worked out exactly. 

The width of a wave-packet describing the spatial motion of a damped harmonic oscillator 
with frequency uoq and mass M in equilibrium with a heatbath is determined by the imaginary 
part of the response function for damped oscillator x"{^) through the FDT, 

= — &) (IS) 

From the discussion by Caldeira and Legget |^, it follows immediatly that this also will 
be the case in equilibrium {eV = 0) for our model, with t] = rjph + rjei. For sufficiently low 
damping, where the energy eigenstates of the oscillator is welldefined, x"{'^) be a narrow 
function centred around ujq, and for low temperature we get. 



(T^ ocl + 2e-^^°/*^^. (14) 

Carrying out a similar analysis in our non-equilibrium case, we find in the same limit, that 

ocl + 2e-'^-«/^^ + (15) 
r]tot{eV)uJo 

Comparing the non-equilibrium expression (p^ with the equilibrium expression (|T^ enable 
us to define an effective temperature by equating the small terms: 

In figure 1 we have plotted the effective temperature as a function of the voltage eV. 

In the case where we can neglect the exponentially small ^-^•^o/kT ^ effective Boltzmann 
factor is proportional to a and thus to eV . The observed rate can now be understood from 
the Arrhenius-factor e"'^"/'^^ dominating the rate-expressions of escape problems in general 
||T6| , where Uq is the barrier- height. If we, as suggested in [^H, assume that f/o ~ S/icuq, we 
conclude that r^^ oc {eVY- In figure 2 the escape rate is plotted as a function of voltage for 
different temperatures. We predict a drastic change in this behavior as the temperature is 
raised which can be verified experimentally. 

We have not, in this calculation, taken into account the dependence of the tunnelma- 
trixelements on the adsorbate position (Tfc(x), Ti{x)). This will contribute to the friction 
as well as the fiuctuating forces on the Xe-atom. Taking this dependence into account also, 
the general scheme described above can be used to calculate the friction coefficient from 
a microscopic starting point, using given expressions for ea(x) and Tfc(x), T;(x). This will 
be applied in a future publication to the problem of laser-induced desorbtion, where the 



resulting temperature dependent friction coefficient is used in the Langevin equation. [|T7 

In conclusion we have outlined a formalism that provide a theoretical foundation of 
the DIMET point of view of desorbtion and applied this to the non-equilibrium Eigler 
experiment . 
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FIGURES 

FIG. 1. The effective temperature as a function of the apphed voltage eV. The absolute ratio 

a/rj is estimated roughly using the ratio, Xe vibrational- to electronic life-time, Tgi/Tph ~ 100, and 
uq S3 3meV. 



FIG. 2. The transfer rate, r , as a function of junction voltage, for T = OK ,4:K ,6K and 8K. 
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